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TECI' "..'.AL MEMORANDUM X-53958 


LARGE ANGLE MnTHOD FOR SPACE VEHIClE ANGULAR 
MOMENTUM DESATURATION USING GRAVITY 
GRADIENT TORQUES 


SUMMARY 


An angular momentum desaturation method is presented for space 
vehicles in orbit. This method utilises the gravity gradient torques and 
therefore avoids the necessity for mass expulsion by reaction jet attitude 
control systems. The desaturation method is based on the assumption that* 
the minimum moment-of-inertia axis has to be aligned with the sun vector 
during orbital daylight for solar experiments. It is further assumed that the 
difference between the two large moments of inertia is negligible. Large 
angular maneuvers of the minimum moment-of-inertia axis during the orbital 
night maximize the gravity gradient torques that are then utilized for 
desaturation of the angular momentum stored by a momentum storage device 
such as control moment gyros. 

The basic idea of the desaturation method is that, for as long as 
possible, the vehicle minimum moment-of-inertia axis keeps a constant 
attitude with respect to a rotating coordinate system in which the local 
vertical is along one of the axes and the orbital angular velocity is along 
another. The constant attitude, with respect to the rotating coordinate system, 
is calculated for propel momentum desaturation. 

To keep a constant attitude in a system rotating with constant angular 
velocity requires, in general, inertial angular rates about all body axes 
simultaneously. These rateo are constant and present no problem for the 
implementation of the desaturation method. 

The gravity gradient torque is perpendicular to both the local vertical 
and the minimum moment-of-inertia axis for the assumed vehicle configuration. 



INTRODUCTION 


Some missions for space vehicles in orbit require an investigation 
of a celestial object which in turn requires an inertially fixed attitude euring 
the time the celestial object is not occiuteu by the earth; e. g. , the Apollo 
Telescope Mount (ATM) and its solar experiments [ 1,2] . In the following, 
ATM nomenclature and definitions will be used; e. g. , sun instead of celestial 
object or night instead of occultation time. etc. The once-considered 
ATM/LM/CSM configuration will also be used, since it had the moment-of- 
inertia distribution typical for the angular momentum desaturation method 
under discussion. 

Gravity gradient, aerodynamic, and other external torques acting on 
the vehicle during the time the vehicle is inertially fixed must be absorbed 
bv an angular momentum storage device, which always has limited capacity. 

A method for momentum desaturation which does not require mass expulsion 
is desired. The fact that an inertially fixed attitude is not required during 
night allows the utilization of the gravity gradient ir. conjvnrtion with vehicle 
attitude maneuvpvo : dCaaturate the momentum storage device. 

Two typically different cases can be identified for the problem of 
momentum desaturation if it is assumed that two of the principal moments 
of inertia of the vehicle are aoout the same and the third is smaller. In the 
first case, the minimum noment-of-inertia axis can be out into the orbital 
plane and kept perpendicular to the sun vector. This case is treated in 
Reference 3 and does not present a problem. In the second case, the minimum 
moment-of-inertia axis is parallel to the sun vector, and a desaturation 
method for this case is treated in Reference 4. However, simulation shoved 
that this method (which worked satisfnetoi iiy lor the gravity gradient torques 
alone) was not effective when significant aerodynamic torques were encountered. 
The method described in the following sections increases the effectiveness by 
maximizing the gravity gradient torques and the time thesp torques are applied. 
This require® the minimum moment-of-inertia axis to h-.ve a constant altitude 
with respect to the gravity gradient direction, and this attitude must be main- 
tained as long ao possible. 



GRAVITY GRADIENT .ORQUE AND MOMENTUM ACCUMULATION 


as 


The gravity gradient torque acting on the vehicle can be expressed 



:>G“ r i r 


(1) 


for a vehicle in circular orbit (which is assumed for further development), 
where 1 is the vehicle iuoment-of-inertia matrix. 12 is the orbital rat?, r 
is a unit vector parallel Ur ilk radius vector from the earth center to the 
vehicle center of mass, and r is defined as 
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When equation (1) is developed for the principal moiuent-of-inertia system, 
P, we have 
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Throughout this report, the assumption is made that there is no 
misalignment between the vehicle principal moment-of-inertia axis and the 
geometric axes. Under this condition the components of the radius vector £ 
with respect to the principal axes are 




where r/^ is the elevation angle of the sun vector from the orbital plane and 
is the orbital angle* from midnight (definitions in Appendix A J. For -~.ero 
altitude error the torque about the vehicle axes becomes 
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If the attitude of the vehicle is held constant between n and n , the 
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following momenta will accumulate (where H is the stored momentum at 
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ANGULAR MOMENTUM DESATURATION METHOD 


The assumed inertia distribution for case two is such that Al and 

x 

Al are almost equal in magnitude, whereas Al is very small. Equation 
y * z 

1 0) shows that only the \-axis has :• bias torque. It is comparable in 



magnitude with the cyclic y -torque. The only way to tie saturate the momentum 
caused by the bias torque is by reversing the sign on the sun elevation, angle 
ij ; i.e. , a large manc.ver about the x-axis at the begUwing and at the end of 

the night is maialatory. Therefore, this maneuver is part of any desaturation 
method for case !«a { I, ij. Visualization of the gravity gradient torque 
(which, for the general case, has Ux& attempted In Reference 6J for case 
two is relatively simple when the large moments of i ne rtia are assumed to be 
equal. The .orp^ s then pc rpendiedar to tae gravity gradient and the 

z -axis (mini mu a moment-of-inertia axis), and it tends to align the latter 

P 

with the former. Then, die gravity gradient desaturation method is developed 
in an iner'iaUy fixed orbital coordinate system. L Abe re the x^-axis is in 

Ifae orbital plane and perpendicular to the sun vector, the ?v-rtx£s is perpen- 
dicular to the orbital plane pointing north, and the z^-axis completes the 
right-handed coordinate system. The z^-axis *oiucjnJes with the projection 
of the sun vector into the orbital plane. 


The x-momentum desaturation is the most critical because, for 
maximum momentum accumulation on the day-half of the orbit (N^l ~ 45 

degrees X inshnb’voiis maneuvers of 90 degrees about the x-axis would be 
required at the terminators to equal the accumulation. Clearly, more capability 
is required; therefore, other maneuvers for x-desaturation are necessary. 

The best that can be done is to keep the z^-axis at a 45-degree angle with 

respect to the gravity gradient direction (maximum gravity- gradient torque# 
and have the projection of the z^-axis into the orbital plane coincide with the 

gravity gradient directi'- m. Under these circumstances the x^-terque is 


proportional to die cosine of the orbital angle, rather than the square 

of the cosine, resulting iu a much larger desaturation capability than the 

maximum day accumulation . The z -axis then describes a cone in inertial 

P 

space. A constant y^-torque results when the z^-axis projection docs not 

coincide with the direction of the gravity g radical but is offset by a constant 
angle. This is used for the y-n> ^centum desaturation. A z-moinenium 
results only when the integration interval is not symmetrical wiUi respect 
to midnight. Xonsymiuetrv is wherefore used for the z- momentum desatura- 
tit n. The desaturation method also has to work when an x-desaturalion (sun 
vector in the oibital plane) is unnecessary and it was found that a tilling of 


G 



the maneuver cone by a small z^- rotation will result in the desired desatura- 
tion capability as shown in the equations developed later. 

A series of single-axis maneuvers is initiated at the beginning of the 
night to arrive at the desired attitude. Constant angular velocities about two 
axes ar? required to hold this constant attitude of the minimum inoment-of- 
inertia axis with respect to the direction of the gravity gradient. Another 
scries of singie-axis maneuvers brings the vehicle back to the desired sun- 
orienled day attitude at the end of the night. The series of maneuvers (always 
about the instantaneous vehicle axes! is given in Table 1. 


TABLE 1. SINGLE -AXIS MANEUVERS 


At Beginning of Night 

| Before End of Night 
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ihe two y-iuaneuvers require further explanation. They are performed in 
the orbital plane (for = pi which is illustrated in Figure 1, The projection 

of the z -axis moves through an angle of 2Ap while the z -axis is coning 
p y P 

between the two sets of single-axis maneuvers. The coning requires an 
inertial angular velocity of the vehicle aliout the y -axis which is resolved 

by u into the y - and z -xves. 
x P P 

All maneuvers are performed by applying precalculated constant 
angular velocities about the appropriate vehicle axes for p recalculated lime 
intervals, identified b\ the equivalent orbital angle »; in Table 2. 










; ' • Figure 1 . Maneuver angles . 


TABLE 2 . ANGULAR VELOCITY COMMANDS FOR 
x, y, AND z VEHICLE AXES 
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with 


K c " ' L/<> L * (11) 

where 12 is the orbital rate and o is the angular rate limit of the vehicle 

JL 

about the large moment-of-inertia axis. 


DESATURATION COMMANDS 


Each desaturation method can be split into two parts. One part is 
concerned with the generation of the desaturation commands; the other part 
is concerned with the scheme to desaturate the commanded momentum. These 
two parts are rather independent of each other; e. g. , several different methods 
for the momentum command generation could use the same method for the 
maneuver angle generation. 

Several facts that affect the momentum command generation must be 
born in mind.: The first is that approximations will have to be made for 
the calculation of the angular momentum desaturated by a given set of 
maneuvers, or the mathematical treatment becomes unmanageable. Secondly, 
the stored angular momentum is the difference between the accumulation and 
the desaturation and does not indicate the desaturation itself. Furthermore, 
the assumptions were made that the large inertias are equal and that there is 
ho misalignment between the principal moment-of-inertia axes and the geo- 
metric vehicle axes. It is therefore impossible to generate an exact command, 
and the desaturation command has to be updated by a summation process such 
that li = x, y, z) 


H Din H Di(n-i) ^ AH Di 


U3) 


where H^. is the present desaturation command, H . , is that of the 

Dm r Dun-1) 

past orbit, and AH is the desired change. The latter has to be chosen 
such that Hp, will home in on the exact solution. Sample data system 
consideration shows the following lorm to be adequate 1 3 . 4 J : 


In 
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n Ain (n-1) Ai(n-l) 


( 13 ) 


i.e. , it is necessary to also consider the past orbit. H^. is an average 
angular momentum dependent on the stored momentum of the following form: 



-7 (H. + H . ) 

2 unax inun 



( 14 ) 


This form averages the peak momenta (during the day-half of the orbit only) 

and provides a selectable offset through H- which allows the conversion, on 

B 

the average, to any point in the momentum space. 

For ease of initialization (initial values for the angles are known, and 
the ones for the momentum commands are not) equation ( 12 ) is replaced by 
equations of the form 



M i(n- 1 ) 


*■ K. AH 
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( 15 ) 


The "memory” is now in the past maneuver commands rather than in the 
summation of the desaturation command changes; cf. equation ( 12 ) . 


DESATURATION MANEUVER EFFECTIVENESS 


The effectiveness of the commanded maneuver angles p » P • P • 

x y z 

Ap . and An has to be established >. th> desired momentum is to be 

y t 

desaturated. The angle Apt . wi*i -.1 — % > be large; p can be large. But 

H r , p_ , and An will be less than i degrees and small angle approxima- 
y z t 

tions can be assumed to simplify the development; i.e. , second-order and 
higher-order terms are neglected. The same reasoning allows neglecting 
any product of the small angles. 


Equation (6) shows the gravity gradient torque about body axes in the 
P- system. The components of v in the P-system are 
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Transformation of the gravity gradient torque to the I-system results in 



Evaluation of equation (lb) with the use of equations (6) and (17) results in 
(second-order terms are neglected; details are in Appendix B) 
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This torque has to be integrated over the orbital angle interval from 

1 ) = A)j - Ap to .» - Ai, + Am . (Tabic 2 and Fig. 1). The aecumu- 

lo l j lb L J 

lated momentum components in the I-system (inertiaily fixed orbital midnight 
system) are (details are in Appendix B) 
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Expressions of the maneuver angles as a function of the commanded 
momentum are needed such that the actual momentum of equations (2u) . (31). 
and (23) equals or approaches the commanded momentum. Additional con- 
straints are needed since there are more maneuver angles than equations. 

The x-moment'un desaturation takes first priority and equation ( 20 ) 
shows that Ap^. should be as large as possible, subject to the available 

maneuver time. It will therefore be made a function of all the other maneuver 
angles and is not considered independent; ef. equation (34). 

Equations (12) and (15) show that there is interest in the gains 
connected with the change c e the maneuver commands. Partial differentiation 
will supply the one for AH : 
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A small positive quantity, k. has to be added to c2p to avoid tho 

possibility for infinite change of p^. which is only en indication of the lack 

ot effectiveness of u . The p -value of the previous orbit can be used for 

x x 

c2p^ since p changes slowly. The \-gain now lias the form 
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l'hc problem has been reduced lo equations ’ml (22; and three 

maneuver angles (p , p . and A/; >. The capability ot p t to desalinate 

a y-momentum is relatively large, and it is ihciMurc logical U* address the 
/-momentum desaturation firs* assuming that ^ can easily compensate tor 
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any effect of p on II , 
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Elimination of p and An from equation (22) with the aid o. equations * 20 ) 
‘l* l 

and (2(>) shows that the actual momen.um, H^. equals the commanded 

momentum, H . Tne relationships of equations (2j) .. <d ( 2 K) also hold 

for the changes of the angles for a chat ge in the momentum command and we 
have 
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With li known, tile actual and the commanded v-momentum can now be* made 
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When Ulc alcove equations arc c ast into ine loc al of equation (15» we have 
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The angular momenta are measured about the vehicle axes (which. under the 
assumption of no misalignment, is the P-syslem; while being sun-oriented, 
l»ut for the development so far, all momenta components were assumed to lie 
in the I-system, and a resolution is necessary: 
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Assume an orbital angle from -p to +p (usually equivalent to the night 
interval; with respect to midnight is available for desaturaticn : then Ap^. is 
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This calculation will be made after the gains K.. are established: i.e. . the 

lj 

gain calculation is performed with the Ap value of the past orbit and Uie new 

Ap is calculated with the maneuver angles of the present orbit; cf. equations (Js) 

through (34j. The reasoning is that the K.^’s are only an approximation. 

whereas a correct Ap is important for maximum cilicier.cy ot the desatura- 
tion method. * 
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CONCLUSION AND RECOMMENDATIONS 


The minimum jicrccntagc of the orbit needed tor full angular momentum 
dcsalurntion us ini; gravity gi ndsr nt torques is :51.9 jxucent, where the following 
assumptions are made: the sunline elevation angle ig. is 4-» degrees (worst case), 
the vehicle has infinite angular rale capability . and no other torques liesuies 
gravity gradient act on the vehicle. The normalized torque profile about the 
xj-axis for this ease is shown n Figure 2 . A finite angular rate causes an 
increase of the minimum desnturation percentage for example, simulation 
shows that the minimum de saturation jH-rcentage ncrcases to about 38 percent 
for a maximum **r.to of 0.01 rad s. Considerate of other external torques 
(aerodynamic, magnetic, etc. ) could further increase the minimum desaiuraiion 
percentage. 


A satellite in a circular orbit with an nllitu-le of 43-1 km. for example. -ias 
a nighttime of only 33.4 percent for 43 dvgr* es. anti the exclusion must lie 
drawn dial valuable daytime which should !ie available for experiments is nec'~>' 
for momentura th "saturation. Unless other means can lie utilized for momentum 
des.it u ration (magnetic torques, etc.), the recommendation is to review the 
mission as to whether it is possible to have the minimum moment -of-inertia 
axis perpendicular to the sunline and simultaneously in the orbital plane (lor 
ease one. see Reference 3). Such an orientation results in very small desaiura- 
lion percentages, osjiecially if it is possible to keep the difference between the 
large moments of inertia small. Angular velocity requirements for case-one ve- 
hicles are about a magnitude smaller than those for ease-two vehicles, allowing 
the use of a smaller capacity for the momentum storage system. 
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APPENDIX A 

DEFINITIONS OF ATM COORDINATE SYSTEMS AND ANGLES 


The angles and coordinate systems used for Uu- AIM aix* shown in 
Figures A-l through A- G. \11 systems were not necessary for the develop- 
ment of the desalu ration method. (The ones not used do not appear in the 
DEFINITION OF SYMBOLS . ) 

T ransf or an lion 


Symbols 


Matrix 

Definition 




The equatorial system is the basic system 
!o which all others are referenced. It is an 

V_ 



earth-centered inertial frame with the 

t 



z-axis toward vernal equinox and the y-axis 

Z E 



toward earth north tFig. A-Hj . 

X u 



The orbital system is an earth-centered 
frame with the z-axis toward the ascending 

v o 

x o 

* *\'*V X E 

node and the y-axis toward orbital north, 
idei. tiffed as a quasi-inertial system 

7 

"O 



(Fig. A-3f. 

X « 



The reference system provides an ideal 
reference for the vehicle. It lias the 

v .« 

X H 

“ '\ ii; , lx o 

z-axis toward the sun center and the x-axis 
in the orbital plane identified as a quasi- 

z ;» 



inc rtial system (Fig. A-5J. 

x v 



l he vehicle system is geometrically 

Y v 

x v 

i,* x K 

ci n to red in the ATM LM CSM with tile 
z-axis toward the ATM. 

S 




* * 
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V 



Tile principal axis system is Mania r lo 
the vehicle system except that the origin 

Y .> 

X i. 

A 


IS at tile* veilleie center oi v.ass and the 
axes arc :*ion„ the principal moment-* *1- 

z i> 



inertia axes. 


1- 



T nuis tor it u lion 
Symbols Matrix 


iXTinition 





X. 



,,; t T ,; y iX 0 


The disturbance system in which aero- 
dynamic and gravity gradient forces are 
defined. Il rotates in Uie orbital plane 
witii tlie vehicle and maintains the z-oxis 
toward earth center. Hie x-axis is 
tangent to the orbit and opposite to the 
velocity vector (circular orbit only) 

(Fig. A-«>. 


Hie sun or ecliptical system is an earth- 
centered inertial frame with the z-axis 
Xg - toward the vernal equinox and the x-axis 

in the plane of the ecliptic (Fig. A -4). 


The intermediate system is an orbital 
coordinate system between tfco O-system 
X - h JX and the R-system. The x-oxis is perpen- 

dicular to Jie sunline and in the orbital 
plane. Th j y-axis is perpendicular to the 
orbital plane pointing north. 


1 !) 
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Figure A -3. ATM orbital coordinate system relative 
to equatorial system X Q - [Xj U y iX fc% . 
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WORTH 



Figure A-4. ATM sun (ecliptic; coordinate system relative 

to equatorial system X' - i4> JX 

z E 





APPENDIX B 

DESATURATION TORQUE AND ANGULAR 
MOMENTUM COMPONENTS 


The components of the desaturation torque should be in an inertial 
frame for ease of integration to obtain the angular desaturation momentum. 
The components of the desaturation torque in the T-system are l see 
equation (li>)J: 
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Evaluation and use of equation (6) leads to 

cij.tr - n r Sft ) - si) U* r r »)♦» r cji 
I I y y x x 't y y x x z x x 

i | cjj.tr - 1 * r sn j - aq.Oi r + r sji ) -r 
z I 't y y x x 't y y x ^x 

'l y y x x . y y x x 


IgT = ^x'z 


Cft 
X X 


or, with equation (17) and neglecting second-order terms. 


T 

-K5T 


-T V 2 Si 


x 


S -M x c2p y s2ij t 

M y U + cij t 

-s2u sn + n c2u (1 - c2n ) 
X t Z X l 


integration of the x-component of T^, from 


„15 - X to , tfi - A, t + X 


results in twith A. - 2 .in A I ) : 

1 x 


23 



» 



”16 


IG 


H Aj = s2fi I c» db) - n u2p | s2i) dr} 
x x J ! t z x ^ t 


7) 


to 


n 


t5 


S2u (sn. - sir _} + p e2u — - 
' '16 'to z x 2 


... ~ e2rj I 
t(> to 


Substitution of n and n and neglecting second-order effects results in 

!-*? tC> 


H 


= s2p SA*X 
x Aj y 


Similarly for the y-eomponent we have 


7 't0 17 to 

H y Aj = p v (f + c 2 #rj | cfc|t s 2 #* x J cij, drjj. 
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to 


»} 


to 


p (1 + e2g ) (ir_ - n ) + p s2p (si) - srr _) 
y x *16 'to z x *t6 'to 


After substitution of the values for tj we have 
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y A| 


M , (1 + c2p J Aju r + ju, s2p _ s A/i 

3 ^ y ^ y 


Development of the z-componenl results in 
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71 to 
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Aj = - S2 u j SI} «i| + a,, c2u | (1 - c2rj J cbj 

Z X * It z \ ~ t I 
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77 to 
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After substitution of the values for rj uc have 
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■A?L s-R sAu + — tf e2u (ii/i - sZlu ) 
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